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Abstract. Truncation-error analysis is a reliable tool in predicting convergence rates of dis- 
cretization errors on regular smooth grids. However, it is often misleading in application to finite- 
volume discretization schemes on irregular (e.g., unstructured) grids. Convergence of truncation 
errors severely degrades on general irregular grids; a design-order convergence can be achieved only 
on grids with a certain degree of geometric regularity. Such degradation of truncation-error conver- 
gence does not necessarily imply a lower-order convergence of discretization errors. In these notes, 
irregular-grid computations demonstrate that the design-order discretization-error convergence can 
be achieved even when truncation errors exhibit a lower-order convergence or, in some cases, do not 
converge at all. 


1. Introduction. These notes are a response to the recently published arti- 
cle [18]. The article applies a truncation-error analysis to evaluate accuracy of finite- 
volume discretization (FVD) schemes on general unstructured grids. The analysis is 
accompanied by computations performed on regular and irregular grids. While we 
agree with the analysis and its conclusions in application to regular smooth grids, 
we consider the truncation-error analysis and some of the conclusions derived from it 
invalid in application to irregular-grid computations. 

On regular grids, convergence of truncation errors (also referred as local errors 
in the literature) is an accurate indicator of convergence of discretization errors (also 
referred as global errors). However, the truncation-error convergence is often mislead- 
ing for FVD schemes defined on irregular (e.g., unstructured) grids. As shown in [18] 
and before this in [17], the second-order convergence of truncation errors for some 
commonly used FVD schemes can be achieved only on grids with a certain degree of 
geometric regularity. Other studies, e.g., [3, 6, 7, 10, 13, 14, 15, 16, 20, 21], showed 
that truncation-error convergence degradation on irregular grids does not necessarily 
imply a degradation of discretization-error convergence. Examples shown in the fol- 
lowing sections confirm that on irregular grids, the design-order discretization-error 
convergence can be achieved even when truncation errors exhibit a lower-order con- 
vergence or, in some cases, do not converge at all. Note that these results do not 
contradict the Lax theorem, which states that consistency (convergence of truncation 
errors) and stability are sufficient (not necessary) for convergence of discretization 
errors. In fact, for some formally inconsistent discretization schemes, it has been 
rigorously proved that the discretization errors converge [3, 6, 7, 13, 16, 21]. 

Article [18] applied a truncation error analysis to FVD schemes for the Poisson 
equation. A thin-layer approximation was analyzed. It was shown that the trun- 
cation error is 0(1) (i.e. , does not converge) in grid refinement unless the grids are 
regular. The discretization error of the scheme was inferred to be non-convergent. By 
coincidence, the particular thin-layer FVD scheme considered in [18] is indeed zeroth 
order accurate on general unstructured grids. The inaccuracy is caused by poor ap- 
proximation to the normal gradients with the thin-layer approximation, not by 0(1) 
convergence of truncation errors. 
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In [18], a general conclusion was drawn that “a compact finite volume approxima- 
tion of the Laplacian has to rely on symmetries in the grid to be first-order accurate.” 
This conclusion is incorrect. For example, a common hnite- volume scheme equiva- 
lent to a Galerkin finite-element approximation on triangles satisfies the definition 
of a compact scheme and is known to have second-order discretization errors (and 
zeroth order truncation errors) on irregular (non-symmetric) grids. We show addi- 
tional computations on general mixed-element grids demonstrating the second-order 
convergence of discretization errors and zeroth order convergence of truncation errors. 

Article [18] also considered a central FVD scheme for an advection equation on 
mixed-element and perturbed quadrilateral grids. Truncation-error analysis showed a 
zeroth-order convergence in the L a c norm. Supporting computations showed a zeroth 
order convergence of discretization errors. It was concluded that FVD schemes for an 
advection equation are non-convergent on non-smooth irregular grids. The conclusion 
is wrong because there are counter examples of FVD schemes with truncation errors 
that do not converge on general irregular grids but with discrete solutions that con- 
verge with at least first order in any norm. The numerical scheme considered in [18] 
is not representative of current practice — the central scheme is known to exhibit 
erratic convergence of discrete solutions because of checkerboard instabilities. Tests 
exposing the erratic convergence of discretization errors of the central scheme are 
shown in Section 4 of these notes. 

In the following Section 2, we provide definitions and discuss relations between 
truncation and discretization errors. Section 3 presents simple one-dinrensional ex- 
amples that expose shortcomings of the truncation-error analysis for irregular-grid 
FVD schemes. Section 4 presents numerical two-dimensional computations for the 
Poisson equation on mixed-element grids and for an advection equation on perturbed 
quadrilateral grids similar to those used in [18]. 

2 . Formulation. For completeness of the presentation, we start from a brief 
formulation of finite-volume discretization (FVD) schemes and definitions of the cor- 
responding truncation and discretization errors. The FVD schemes are derived from 
the integral form of a steady conservation law 

/(F-fi) dT = JJfdn, (2.1) 

r n 

where / is a forcing function, D is a control volume with boundary F, n is the outward 
unit normal vector, and F is the flux vector that may include viscous and/or inviscid 
contributions. The general FVD approach requires partitioning the domain into a set 
of non-overlapping control volumes and implementing numerically equation (2.1) over 
each control volume. In this paper, we are considering node-centered FVD schemes, 
where solution values are stored at the mesh nodes. 

For two-dimensional (2D) node-centered FVD schemes, a median-dual partition 
is constructed by connecting the mass centers of the primal-mesh cells with the mid- 
points of the surrounding edges (Figure 2.1). These non-overlapping control volumes 
cover the entire computational domain and compose a mesh that is dual to the primal 
mesh. The discretization is applied at a sequence of refined grids composed of trian- 
gular and quadrilateral cells satisfying the consistent refinement property [8, 20]. The 
property requires the maximum distance across the cells to decrease consistently with 
increase of the total number of grid points, N. In particular, the maximum distance 
should tend to zero as V -1 / 2 in 2D computations. 
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Fig. 2.1. Median- dual partition for node- centered finite- volume discretizations. Pq — P 4 denote 
grid nodes. The control volume (dual cell) around Pq is shaded. 


The main accuracy measure of any FVD scheme is the discretization error , Ed , 
defined as the difference between the exact discrete solution, Q h , of the discretized 
equations (2.1) and the exact continuous solution, Q, to the differential conservation 
law 


VF = /, (2.2) 

E d = Q — Q h \ (2.3) 


Q is sampled at grid nodes. 

Truncation error, E tl measures the accuracy of the discrete approximation to 
the differential equations (2.2). For FVD schemes, the traditional truncation error 
is usually defined from the time-dependent standpoint [19, 22]. In the steady-state 
limit, it is defined (e.g., in [10]) as the residual computed after substituting Q into 
the normalized discrete equations (2.1), 


E t = 


1 

W\ 



(2.4) 


where F h is a reconstruction of the flux F at the boundary T, f 1 is the measure of 
the control volume, 


|fi| = JJ dSl, (2.5) 

n 

f h is an approximation of the forcing function / on fl, and the integrals are computed 
according to some quadrature formulas. 

For a given FVD scheme, we define a residual function 


R(q) = 


1 

W\ 



( 2 . 6 ) 
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where q is a discrete function defined at the grid nodes. The residual accounts for 
interior discretization and all boundary conditions. Note that by definition 

R(Q h ) = 0 (2.7) 


and 


R(Q) = E t . (2.8) 

Substituting (2.3) into (2.7), 

R{Q-E d )= 0, (2.9) 


and assuming the discretization error to be small comparing to the exact continuous 
solution, Q, ( \Ed\ « |Q|), we obtain 


E t - J(Q)E d « 0, 


( 2 . 10 ) 


where J (Q) is the Jacobian of R(q) computed for q = Q, 


J (Q) 


J_d_ 
|f2| dq 


/ 


.r 


(F h (Q) • n) dT . 


Thus, the discretization error can be estimated as 


E d « J~ X E U 


( 2 . 11 ) 


( 2 . 12 ) 


where the inverse Jacobian is a convolution operator with a kernel that usually pos- 
sesses certain smoothness properties. The discretization error in a point can be ap- 
proximated as a weighted sum of the undivided truncation errors, 

E d KY / w j (E t ) j \il j \, (2.13) 

3 

where summation index j goes over all control volumes and weighting factors, Wj , are 
determined by the convolution kernel. 

Two norms are used in this paper to evaluate magnitudes of the discrete errors: 
the Loo norm (also referred as max-norm) is computed as the maximum absolute 
value of a discrete function; the L\ norm is computed as the mean absolute value, 
i.e. , the sum of the absolute values divided by the number of discrete entries. Thus, 
either norm of a constant discrete function is the same constant. 

Irregular grid computations, in which convergence of discretization errors sur- 
passes the convergence of truncation errors, were reported and analyzed in literature 
since the 1960s (see, e.g., [3, 6, 7, 10, 13, 14, 15, 20, 21]). In [13], discretization schemes 
demonstrating such convergence were called supra- convergent. While a rigorous proof 
of discretization error convergence for FVD schemes on general irregular grids is not 
yet available, there are several recent publications addressing supra-convergence on 
irregular grids. Papers [7, 16] consider formally inconsistent (no truncation-error con- 
vergence) schemes for advection equations in ID and 2D and prove convergence of 
discretization errors; Barbeiro [3] proves second-order convergence of discretization 
errors for “inconsistent” discretizations of 2D elliptic equations on nonuniform grids. 
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A qualitative explanation of supra-convergence observed for FVD schemes on 
irregular grids was suggested by Giles [10]. Giles’ arguments can be illustrated as 
follows. First consider a node-centered FVD scheme on a sequence of regular grids. 
On such grids, truncation and discretization errors would converge with the same 
design order. Now consider effects of a local irregularity introduced, for example, by 
shifting one node. The node shift leads to changes in faces (and corresponding fluxes) 
of the control volume containing this node and its neighboring control volumes. Fluxes 
outside of this immediate neighborhood remain unchanged. Each (changed) flux at 
an internal face contributes equal and opposite amounts to the undivided truncation 
errors on both sides of the face. Thus, the change introduced into the local truncation 
error can be large (the change in a flux divided by the control volume), but the mean 
value (and, often, higher moments) of the changes introduced into the undivided 
truncation errors is zero. If the convolution kernel is sufficiently smooth, e.g., for all 
nodes outside of the immediate neighborhood of the perturbed node, the weighting 
factors do not vary too much over the control volumes with perturbed truncation 
errors, then the discretization error changes are close to zero outside of this immediate 
neighborhood. Even within this neighborhood, the discretization error changes are 
small comparing with the changes in the truncation errors because of the averaging 
effect. 

3. One-dimensional examples. Examples in this section show that common 
FVD schemes on one-dimensional (ID) irregular grids exhibit supra-convergence. Sim- 
ilar computations have been demonstrated in [15, 16]. 

3.1. Spatial discretization grids. A ID discretization grid is defined as a com- 
bination of the primal and dual nodes. The solution values are located at the primal 
nodes; the fluxes are located at the dual nodes. For node-centered discretizations, 
a natural strategy is to place the primal mesh first and, then, use this mesh as a 
reference for placing dual control volumes. 

The discretization grids used in this section are designed to demonstrate effects 
of grid irregularities and are described as follows. The first and the last of the N + 1 
primal nodes, Xi, i = 0, 1, . . . , N, are always located at the ends of the computational 
interval; the interior nodes can be distributed either uniformly or randomly. Either 
distribution retains the nodal ordering and ensures that the maximal distance between 
the neighboring nodes is 0(1/N). Let s,, i = 0, 1, . . . , (N+ 1) denote the flux locations. 
The first and the last fluxes are also located at the ends of the interval. The location 
of an interior flux, Sj, is always between the primal nodes Xi-i and aq and initially 
defined as the midpoint; the dual node ,s,; may then be randomly perturbed about the 
midpoint. 

Specifically, on an interval x £ [a, b ], the primal nodes are distributed according 
to 

x 0 = a; Xi = a + (i + > * = 1, . . . , (N - 1); x N = b; (3.1) 

where r, is either zero (uniform primal mesh) or a random number —0.4 < r, < 0.4 

(random primal mesh). The dual-mesh nodes are computed accordingly as 

s 0 = a; Si = x^ + di(xi - Xi-i) , i = l,...,N; s N+ i=b; (3.2) 

where di = 0.5 and di = 0.5 + rf correspond to unperturbed and perturbed dual 

meshes, respectively; here rf is a random number —0.25 < rf < 0.25. Grid examples 
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(a) Uniform primal mesh; unperturbed dual mesh. 


* 


(b) Uniform primal mesh; perturbed dual mesh. 


• i •— i •— i • i • I • i— • 


(c) Random primal mesh; unperturbed dual mesh. 


Fig. 3.1. Examples of one- dimensional discretization grids: black bullets denote primal mesh 
nodes, vertical tic-marks denote dual mesh nodes. 


are shown in Figure 3.1. The computations in this section refer to tests performed 
on the interval x 6 [0, 1] using a sequence of grids with the total number of primal 
grid nodes increasing as N = 2 n + 1; n = 3, 4, . . . , 14. All random grids are generated 
independently, so there is no regularity in the limit of grid refinement. The equivalent 
mesh size is taken as h = -k. 



Equivalent Mesh Size 


Fig. 3.2. Convergence of discretization and truncation errors in grid-refinement computations 
for the constant- coefficient convection equation. The tests are performed with random primal meshes 
and unperturbed dual meshes. 
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3.2. Steady convection equation. The first example is a steady constant- 
coefficient convection equation 

d x U = f(x), U( 0)=ip, (3.3) 

which is satisfied with the exact solution U = sin (x) , f = cos(x). FVD equations 
are formed as follows 


F Si+1 ^ F Si = sin(s i+ i) - sin(sj), i = 0, ...,7V. (3.4) 

The fluxes, F s . , approximate solution values at the flux locations, .sy, and are computed 
by fully-upwind extrapolations from the primal nodes (except for the first interior dual 
node) as 


F s o = ip, 


_ (si—xq)Ui + (xi — si)Uq 

K, = 


F S1 = 


Xl—Xo ’ 

(si-Xi- 2 )Ui-l-(si-Xi-l)Ui-2 


i = 2,...,N + l 


(3.5) 


where [/,; is a discrete approximation to U(xi). Boundary conditions are implemented 
weakly by defining the inflow flux. The forcing term is integrated exactly over the 
dual volume. 

Figure 3.2 shows convergence rates obtained in grid-refinement on irregular dis- 
cretization grids involving random primal meshes and unperturbed dual meshes. In 
either the or L\ norm, truncation and discretization errors converge with different 
orders. 


Primal 

Dual 

Grid-refinement computations 

Mesh 

Mesh 

Discr. Error 

Trunc. Error 

Uniform 

Unperturbed 

~o[ F) 

0(h z ) 

Uniform 

Perturbed 

om 

0(h) 

Random 

Unperturbed 

0(W) 

0(h) 

Random 

Perturbed 

om 

0(h) 


Table 3.1 

Convergence of discretization and truncation errors for the constant- coefficient convection equa- 
tion. 


Table 3.1 summarizes discretization and truncation error convergence rates ob- 
served in computations. The results demonstrate that grid irregularity strongly affects 
the truncation-error convergence, but has no effect on the convergence order of the 
discretization errors. 

3.3. Time dependent convection equation. The time-dependent constant- 
coefficient convection equation is defined on the domain (t, x) € [0, 1] x [0, 1] as 

d t U + d x U = f{x), U(t, 0) = ip(t), U(0,x) = <p(x), (3.6) 

and is satisfied with the exact solution U(t,x) = sin (2n(x — t)) , / = 0. In com- 

putations, the number of time levels is equal to the number of spatial grid nodes, 
N + 1. The time step is constant, ht = 1 /TV. The same spatial grids are used at all 
time levels. 

FVD equations at time level k are formed as follows 

(-/f + (s i+ i - 8i ) + Ft +i - F* = 0, i = 0, . . . , TV. 
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(3.7) 




Here, is a discrete time derivative and is a discrete approximation to the 

forcing term; both are defined at grid node X{. For k > 1, is the second-order 

backward difference scheme, 


d t u i = ^ - 2t/T 1 + , (3.8) 

and, for k = 1, the time derivative is approximated by the first-order backward 
difference scheme, 


d h t U\ = ^(Ut- U?) . (3.9) 

The spatial fluxes, F s fc . , are formed according to (3.5). The time dependent truncation 
error is evaluated for the exact solution, U, at each dual volume as 

jpk j?k 

E t = -f* + d?U+-^ -. (3.10) 

Sj+l — Si 



Fig. 3.3. Convergence of discretization and truncation errors in grid-refinement computations 
for the time- dependent constant- coefficient convection equation. The errors are computed at the 
final time t = 1. The tests are performed with random primal meshes and unperturbed dual meshes. 

Figure 3.3 shows convergence rates obtained in successive grid-refinement, which 
doubles the number of primal mesh intervals and time steps. Spatial grids involve 
random primal meshes and unperturbed dual meshes. The coarsest mesh includes 9 
primal nodes (8 intervals) and uses 9 time levels (8 steps). Truncation and discretiza- 
tion errors are computed at the final time t = 1. In either the L 0 0 or L\ norm, the 
truncation errors converge with the first order and discretization errors converge with 
the second order. 


3.4. Diffusion equation. The third set of one-dimensional tests is performed 
for the constant-coefficient diffusion equation 

d xx U = f(x), U( 0 )=^ 0 , U{ 1 ) = - 01 , ( 3 - 11 ) 

which is defined on the interval x £ [0, 1], with the exact solution U = sin(a:), 
/ = — sin(a’). FVD equations are formed at the interior dual volumes as 

Fs i+ i — Fsi = cos(sj+i) — cos(sj), i = (3-12) 

Fluxes approximating the solution derivative are defined as 

F Si = Ui ~ Ui ~\ i=l,...,N, (3.13) 

•%i %i— 1 

where Ut is a discrete approximation to U(xi) and the forcing term is integrated 
exactly. Dirichlet boundary conditions are strongly enforced as 

U 0 = ipo, U N+1 = </h. (3.14) 


Primal 

Dual 

Grid-refinement computations 

Mesh 

Mesh 

Discr. Error 

Trunc. Error 

Uniform 

Unperturbed 

om 

0(F) 

Uniform 

Perturbed 

0(F) 

0(1) 

Random 

Unperturbed 

O(F) 

o(h) 

Random 

Perturbed 

O(F) 

0(1) 


Table 3.2 

Convergence of discretization and truncation errors for the diffusion equation. 


Table 3.2 summarizes the discretization and truncation error convergence. The 
main observation is that the discretization errors converge with second order even for 
those irregular grids for which truncation errors are 0(1). 

4. Two-dimensional computations. In this section, we demonstrate 2D com- 
putations for the Poisson equation and for a constant-coefficient advection equation 
on general irregular grids. More details and other computations, including inviscid 
flow equations on general irregular grids, can be found in [1, 8, 20]. 

4.1. Poisson equation. A widely used second-order FVD scheme is shown here 
for the Poisson equation, 


A U = f, 


(4.1) 


with Dirichlet boundary conditions. In the FVD scheme, the conservation law 


/ 


VU-ndr = 



(4.2) 


is enforced on node-centered control volumes constructed by the median-dual parti- 
tion; n is the outward unit normal to the control-volume surface. In distinction from 
the thin-layer approximation considered in [18], the scheme uses the Green-Gauss 
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Fig. 4.1. Illustration of gradient reconstruction for viscous terms on mixed grids with median- 
dual partition. 


approach to provide accurate approximations to the normal gradients at the control- 
volume boundaries. A series of refined mixed-element grids composed of triangles and 
quadrangles is considered. 

With reference to Figure 4.1, the integral flux through the dual faces adjacent to 
the edge [Po,Pi] is computed as 


/ VU-ndr = VU fl -n R + VU L n L , (4.3) 

ABC 

where iir and riL are directed areas of the corresponding dual faces. The gradient is 
reconstructed separately at each dual face as follows. 

For the triangular element contribution, the gradient is determined from a Green- 
Gauss evaluation at the primal-grid element, 

VU l =VU 0 i 4 . (4.4) 

The gradient overbar denotes a gradient evaluated by the Green-Gauss formula on 
the primal cell identified by the point subscripts. With fully-triangular elements, 
the formulation is equivalent to a Galerkin finite element scheme with a linear basis 
function [2, 4]. For the quadrilateral element contribution, the gradient is evaluated 
as 

VU fl = VU0234 + [ Ui ~ U °. - VU0234 • e 04 l e 04 , (4.5) 

L l r 4 — r o| J 

where r , is the coordinate vector of node P,; and 


eo4 


i~4 - r 0 

| r 4 - r 0 | 


(4.6) 


is the unit vector aligned with the edge [Pq, P 4 ]. Note that for grids with dual faces 
perpendicular to the edges, the edge-gradient is the only contributor. This approach 
to the gradient reconstruction is used to decrease the scheme susceptibility to odd- 
even decoupling [11, 12]. In all cases, the linear solution reconstruction leads to a 
first-order flux (gradient) reconstruction accuracy. 
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(a) A representative 17x17 mixed-element 
unstructured grid. 



(b) Convergence of the discretization and 
truncation errors. 


Fig. 4.2. Convergence of the discretization and truncation errors for the Poisson equation 
solved on mixed- element unstructured grids. 


A sequence of globally refined mixed- element grids is generated on the unit square 
(x, y ) € [0, 1] x [0, 1] with 9, 17, 33, 65, and 129 points in the horizontal and vertical 
directions. Each grid is formed from an underlying Cartesian grid. Grid irregularity is 
enforced independently on each grid by random splitting (or not splitting) Cartesian 
cells into triangles; approximately half of the quadrangles are split. Figure 4.2(a) 
shows a representative 17 x 17 grid. 

The exact solution is taken as U = sin(7rx) + sin(7n/), and the corresponding forc- 
ing term is / = — 7r 2 [sin(7r;r)-|-sin(7ri/)]. The L\ norms of truncation and discretization 
errors are shown in Figure 4.2(b) versus an effective mesh size, taken as the L\ norm 
of the square root of the dual volume. Discretization errors converge with the design 
second order; truncation errors do not converge. Although not shown, the error con- 
vergence rates in the L ^ norm are the same as the Li-norm rates. On regular (either 
triangular or quadrilateral) meshes, both the truncation errors and the discretization 
errors converge with second order. In [18], 0(1) convergence of truncation errors on 
irregular grids was interpreted as an indication that the corresponding discrete solu- 
tions do not approximate the continuous ones; this example clearly shows that such 
inference can be wrong. The 0(1) convergence of the thin-layer scheme discretization 
errors shown in [18] is caused by poor approximation to the normal gradients provided 
by edge-derivatives on skewed grids, not by 0(1) convergence of truncation errors. 

4.2. Advection equation. A constant coefficient 2D advection equation 


(a • V)f7 = /, (4.7) 

is defined on the unit square (x, y) G [0, 1] x [0,1]; a = (a, b) T — (1, 1) T is a constant 
vector. The characteristic direction is defined as £ = (ax + by)/c where c 2 = a 2 + b 2 . 
The cross-characteristic direction is defined as r] = (— bx + ay)/c. 

There are three exact solutions considered: 


U\ = sin(47rc£); f± = 47rc 2 cos(47tc£). 
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(4.8) 


u 2 = c (£ + £ 2 / 2 ); 


f 2 = c 2 (l + 0- 


(4.9) 


U 3 = sin(7rc?7); f 3 = 0. 


(4.10) 


The first two solutions have variation in the characteristic direction; the third solution 
has variation in the cross-characteristic direction only. 

In the edge-reconstruction FVD schemes considered in this paper, the conserva- 
tion law 


j) (a • h)U dT 



(4.11) 


is again enforced on node-centered control volumes constructed by the median-dual 
partition. The forcing term integration over the control volume is approximated as 
the node value multiplied by the volume f!| . The term edge-reconstruction empha- 
sizes that the quadrature scheme used for computing the integrals over the dual-cell 
boundaries employs fluxes reconstructed at the midpoints of the edges connecting the 
grid nodes. The numerical flux through the faces adjacent to the edge [P 0 ,Pi] is 
computed as 

F h = (a • n)U, (4.12) 

where U is a solution reconstruction at the edge midpoint and n is the combined- 
directed-area vector n = nL + n^, where nL and iir are outward normal (directed- 
area) vectors of the left and right segments, respectively; see Figure 4.3 for illustration. 
FVD schemes with first- and second-order upwind reconstructions, which are scalar 
counterparts to widely used Euler schemes [1, 2, 4, 5, 10, 11, 12, 14, 17, 20], are 
considered as well as the central scheme used in [18]. 



Fig. 4.3. Illustration for edge-reconstruction flux integration scheme in the interior. 

The first-order upwind edge-reconstruction scheme approximates U from the so- 
lution at the upwind node. The second-order accurate upwind-biased (Fromm-type) 
edge-reconstruction scheme uses the least-square gradient at a node to reconstruct 
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the solution at the edge midpoints. Second-order convergence of discretization errors 
has been demonstrated on general (e.g., irregular non-smooth) triangular grids. On 
irregular quadrilateral and mixed-element grids, the scheme accuracy deteriorates to 
first order; examples and explanations can be found in [8, 20]. 

Boundary conditions are enforced weakly through the boundary fluxes. At the 
inflow boundary, U is specified from the exact solution. At the outflow boundary, 
the boundary flux is determined from the solution value at the nearest node point, 
referred to here as a weak closure condition. This closure is approximately a first-order 
upwind differencing closure and is also that advocated in [18]. For triangular cells at 
the boundary, a special flux integrations scheme is applied to provide second-order 
accuracy; details can be found in [2, 8, 20]. 

4.2.1. Perturbed quadrilateral grids. This section reports on investigation 
of the claim made in [18] that edge-reconstruction FVD schemes for perturbed quadri- 
laterals are zeroth order. 

The considered uniform quadrilateral grids are Cartesian grids (Figure 4.4 (a)) 
defined on the domain x G [0,1] and y G [0,1]. For the perturbed quadrilateral 
grids, all grid points are perturbed in both x and y directions by random shifts in the 
range [—0.1, +0.1] of the local mesh size (Figure 4.4 (b)). The sequences of globally- 
refined grids are typically generated with 2 n + 1 points in both directions, where 
n = 3, 4, 5, 6, 7, 8. 



1 

0.8 

0.6 

0.4 

0.2 

0 



X 


(a) Uniform grid. 


(b) Perturbed grid. 


Fig. 4.4. Typical computational grid (shown with 24 points in each direction of the underlying 
Cartesian mesh). 

In the figures that follow, the L\ norms of truncation and discretization errors are 
shown versus an effective mesh size parameter, taken as the L\ norm of the square root 
of the dual volume. The filled symbols denote truncation errors and the open symbols 
denote discretization errors. Circles, squares, and triangles refer to simulations with 
the first, second, and third exact solutions, respectively. 

The results for the first-order accurate upwind edge-reconstruction scheme are 
shown in Figure 4.5. For the uniform grids, both truncation errors and discretization 
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(a) Uniform grid. 


(b) Perturbed grid. 


Fig. 4.5. Errors in grid refinement for first- order upwind edge- reconstruction scheme. Circles 
correspond to solution U\, squares correspond to solution U2, and triangles correspond to solution 
U3. Dashed line indicates the first-order slope. 


errors show a first-order convergence. For the perturbed grids, the truncation errors 
show one-order reduction in convergence, but the discretization errors remain first- 
order accurate. Variations of the magnitude of the discretization error norm between 
the uniform and perturbed grids with the same number of points are very small. 

Figure 4.6 shows the results for the second-order accurate upwind-biased (Fronnn- 
type) edge-reconstruction scheme. Again, for the uniform grids, both truncation errors 
and discretization errors show a second-order convergence. For the perturbed grids, 
the truncation errors show no convergence, i.e., zeroth-order convergence. On coarse 
grids, the discretization error convergence is better than first order. For solutions U 2 
and U 3 , the discretization errors asymptote clearly to first-order accuracy on finer 
grids. The results shown in Figures 4.5 and 4.6 clearly contradict to the conclu- 
sion drawn in [18] that FVD schemes that exhibit the zeroth order convergence of 
truncation errors are not convergent on irregular grids. 

For illustration of erratic behavior of central schemes, the results for the second- 
order accurate central scheme used in [18] are shown in Figure 4.7; this discretization 
reconstructs the solution by averaging the two nodal values on either side of an edge. 
This scheme is susceptible to checkerboard instabilities and boundary conditions are 
prominent in eliminating unstable modes. For the uniform grids, both truncation 
errors and discretization errors show a second-order convergence. Although not shown, 
the discretization errors exhibit a checkerboard pattern over the mesh. 

For the perturbed grids, the truncation errors show no convergence, although we 
do not show their variation. The variation of the discretization errors is shown for two 
grid sequences. The sequences are derived from 9x9 and 13 x 13 coarsest grids. The 
discretization errors are quite erratic with grid refinement, as is typical of schemes 
with checkerboard instabilities. It might be inferred from the few coarsest grids that 
the discretization errors associated with U 2 and U 3 do not converge, but continued 
refinement shows that the discretization errors are converging, but apparently at a 
rate slower than first-order. 
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(a) Uniform grid. 


(b) Perturbed grid. 


Fig. 4.6. Errors in grid refinement for second-order Fromm-type edge-reconstruction scheme. 
Circles correspond to solution U\, squares correspond to solution U2, and triangles correspond to 
solution U3 . Dashed and dotted lines indicate the first and second order slopes, respectively. 




(a) Uniform grid. 


(b) Perturbed grid. 


Fig. 4.7. Errors in grid refinement for second-order central edge-reconstruction scheme. Circles 
correspond to solution Ui, squares correspond to solution U2, and triangles correspond to solution 
U3. Dashed and dotted lines indicate the first and second order slopes, respectively. 


Numerical results presented in this section show both the fully-upwind first-order 
edge-reconstruction scheme and an upwind-biased Fromm-type edge-reconstruction 
scheme to have first-order convergence of discretization errors on perturbed quadri- 
lateral grids. These results are consistent with the results reported in [1, 8, 20] for 
inviscid flow equations. For the central edge-reconstruction scheme, the discretiza- 
tion errors are highly oscillatory spatially and the convergence of discretization error 
norms is highly erratic. The convergence of the discretization errors is apparently 
less than first order, but nonetheless the errors converge in all norms. The erratic 
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behavior of the discretization errors and insufficient grid refinement may have led to 
an incorrect conclusion in [18]. 

Remark on mixed grids computations. Finite-volume computations on general 
mixed-element grids have been demonstrated in [9, 20] for the Euler equations. The 
results showed 0(1) convergence of truncation errors and first order convergence of 
discretization errors in all norms. On mixed-element grids, with non-smooth interfaces 
between two regions of regular grids, such as in [18], convergence in integral norms 
was even better, approaching second order in the L 2 norm. Zeroth order convergence 
of the Loo norms of the discretization errors and partial-order convergence of the 
L 2 norms of the discretization errors reported in [18] for advection computations on 
mixed grids is attributed to insufficient refinement and the erratic convergence of 
discretization errors characteristic for central schemes with checkerboard instabilities. 

5. Conclusions. Truncation-error analysis is a reliable tool in predicting con- 
vergence rates of discretization errors on regular smooth grids. However, it is often 
misleading in application to hnite- volume discretization schemes on irregular (e.g., 
unstructured) grids. Convergence of truncation errors severely degrades on general 
irregular grids; a design-order convergence can be achieved only on grids with a cer- 
tain degree of geometric regularity. Such degradation of truncation-error convergence 
does not necessarily imply a lower-order convergence of discretization errors. In these 
notes, irregular-grid computations with finite- volume discretization schemes demon- 
strate that the design-order discretization-error convergence can be achieved even 
when truncation errors exhibit a lower-order convergence or, in some cases, do not 
converge at all. 
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